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Abstract 

Schwartzenberg et al. 1995| , jlmpey et al. 1996| , pal 
We ha ^e used a large lormat UUP camera to survey canton et al. 1997| , |Kashikawa et al. 1998|) . In this 



the nearby Fornax cluster and its immediate environ- 
ment for low luminosity low surface brightness galax- 
ies. Recent observations indicate that these are the 



paper we extend this survey work to a study of the 
nearby Fornax cluster and its immediate environment 



most dark matter dominated galaxies known and so 
they are likely to be a good tracer of the dark matter 
in clusters. We have identified large numbers of these 
galaxies consistent with a steep faint end slope of the 
luminosity function (a ~ —2) down to Mb ~ —12. 
These galaxies contribute almost the same amount 
to the total cluster light as the brighter galaxies and 
they have a spatial extent that is some four times 
larger. They satisfy two of the important predictions 
of N-body hierarchical simulations of structure for- 
mation using dark halos. The luminosity (mass ?) 
function is steep and the mass distribution is more ex- 
tended than that defined by the brighter galaxies. We 
also fii|id a laige cuiicenlialiuii of luw suiface biighL- 



(see also Caldwell 1983, Davies et al. 1988, Fergu- 



son 1990, Irwin et al. 1990, Bothun et al. 1991 



Held fc Mould 1994, Morshidi-Esslinger et al. 1999a 



X 



a laige cuiicenlialiuii ui luw sui 

ness galaxies around the nearby galaxy NGC1291. 

Keywords: cosmology: dark matter: galaxies: indi- 
vidual (NGC1291), dwarf galaxies, luminosity function: 
clusters individual (Fornax) 



Morshidi-Esslinger et al. 1999b ). 

Two fundamental cosmological issues are the spa- 
tial scale over which mass is distributed in the Uni- 
verse and the mass spectrum of the collapsed objects 
within the Universe. Given that dark matter dom- 
inates the mass the challenge is to measure these 
quantities for the dark, not visible, matter. One 
obvious solution is to measure the spatial distribu- 
tion and mass spectrum of the most dark matter 
dominated galaxies we know of: these are the LSB 
and dwarf galaxies ( McGaugh fc de Blok 1998|, Irwin 



fc Hatzidimitriou 1995, Carignan fc Freeman 1988 



Dckci & ssnnggD . 



Historically the intrinsically bright and high sur- 
face brightness galaxies have been used to define the 
spatial scale and distribution of structures in the Uni- 
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verse (see for example Maddox et al. 1990 , ia Costa 
et al. 1998j ). Their clustering length scale, as mea- 



Over the last few decades new populations of galax- 
ies have revealed themselves through the use of high 
sensitivity detectors at dark observing sites. From 
dwarfs to giants, galaxies of Low Surface Brightness 
(LSB) are a strikingly new addition to the galactic 
inventory of the Universe ( Davies et al. 1998 , Impcy| 
|fc Bothun 1997|). Having identified this new class of 



object we have now moved on into the era of defin- 
ing their cosmological importance in relation to the 
brighter, more well known galaxies. To do this it 
has been necessary to survey significant areas of sky 



sured by such techniques as the two point spatial 
correlation function, defines the bright galaxy mass 
scale. The observed clustering scale is smaller than 
that inferred from simulations of structure growth 
using currently popular CDM models of dark haloes 
(see for example [Frenk et al. 1996 , Jenkins et al. 1998) 
and so it is often assumed that the bright galaxies are 
biased tracers of the underlying mass ( Pekel fc Rees| 
Frenk et al. 199(f ). In these models the bright 
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at LSB levels (e.g. [Morshidi-Esslinger et al. 1999a , 



galaxies are more concentrated than the dark matter 
because they form at the peaks of the initial density 
fluctuations. They also undergo more mergers than 
the objects in the lower density outer regions of clus- 
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ters and so they can become more 'visible' due to 
enhanced star formation. In this paper we intend to 
show that, for one nearby cluster at least, there is a 
population of very LSB low luminosity galaxies that 
contribute about the same luminosity to the cluster 
as the brighter galaxies, but which defines a much 
more extended spatial scale than the bright galaxies. 
Our expectation is that these galaxies are dark mat- 



ter dominated (cf. the "dark galaxy" of Carignan & 
Freem; m 1988[) , having mass-to- light ratios of 100 or 
more similar to the mass-to-light ratios determined 
for galaxy clusters ( David et al. 199*51 ) . 

It is also difficult to relate the mass spectrum in- 



ferred From numerical simulations of structure growth 

(for example CDM) to what we observe. The sim- 
ulations model the mass spectrum of the dominant 
dark matter haloes yet what we observe is the lu- 
minosity function of galaxies. These two are related 
through the total mass-to-light ratio of each object, 
but the mass-to-light ratio is almost certainly a func- 
tion of galaxy luminosity, type etc. and so relating 
luminosity to total mass is not easy. In addition, 
various forms of evolution (mergers, tidal stripping, 
star formation, galactic winds etc.) may have al- 
tered the mass spectrum over time. To overcome 
these problems numerous N-body simulations incor- 
porating ever more complicated physics have been 
carried out in an attempt to relate the initial condi- 
tions ( Press Schector 1974| ) to what we observe (see 
for example [Frenk et al. 1996] , |Kauffman et al. 1997] ) . 
Even with these complications the models of hierar- 
chical structure formation with galaxy evolution still 
predict a steep faint end slope to the local mass func- 
tion, similar to the presumed initial conditions. This 
appears to be a rather robust conclusion. The local 
cluster galaxy mass function predicted by the models 
has a low mass slope of w -2 , flFrenk et al. jggp . 
This is much steeper than that generally found for the 
luminosity function of field galaxies ( Loveday et al.| 
199 2|, |Marke et al. 1994] but see also |Loveday 199^ 
However there are a number of recent determinations 
of a steep luminosity function for clusters ( |Phillipp"s| 



{M B « -12). 

Our prime objective in this paper is to compare 
the spatial distribution and total luminosity of three 
samples of Fornax cluster galaxies selected by their 
surface brightness. We define a bright galaxy sam- 
ple typical of that used to measure the two point 
spatial correlation function and the galaxy luminos- 
ity function, a LSB sample typical of recent work to 
determine the numbers of LSB dwarf galaxies in clus- 
ters and a very LSB (VLSB) sample typical of what 
can now be detected using a large area CCD. The 
bright and LSB samples have been discussed previ- 
ously by us ( pisney et al. 1990| , |Morshidi-Esslingcr 
et al. 1999a] b, henceforth MDS). The VLSB sample 



et al. pJ98|, |Valotto et al. 1997| , |Smith et alT 1996| , 
Bernstein ct al. 1995T ). We will find observationally 



for the Fornax cluster a steep cluster luminosity func- 
tion consistent with the expected mass spectrum in- 
ferred from N-body simulations and show that the 
steep slope extends to intrinsically faint magnitudes 



is new and so we will start by describing these data. 

2 Observation and image process- 
ing of the VLSB sample 

The data were obtained using the Curtis Schmidt 
96cm Telescope at the CTIO during the nights of 
21-24th of October 1997, using the 2048x2048 CCD 
with 2.3 arc sec pixels. The total area observed 
was 13.8 sq deg, forming a strip ~9 deg long out- 
wards from the cluster center (NGC1399) towards 
the lenticular barred galaxy NGC1291 (figure [l]). 
We chose this particular strip of sky because MDS 
(the LSB sample, see below) had identified numerous 
quite large LSB galaxies in this area of sky around 
NGC1291. 

There are 13 individual fields in the strip, each of 
1.3° x 1.3°. Each field consists of four exposures of 
1050sec, median stacked to give good signal to noise, 
while keeping saturation to a minimum. The different 
exposures were "jittered" around the nominal field 
center to improve flat fielding. 

2.1 Data reduction 

The CCD uses two read out amplifiers, so all steps of 
data reduction were performed twice, once for each 
amplifier. The original images (data, bias frames, 
flats and standards) were cut into two parts (accord- 
ing to the reading amplifier) and then each side was 
processed separately. The data remained separated 
for the whole data reduction and detection process. 

Prior to median stacking, the image data were bias 
subtracted and flat fielded using the standard tech- 
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Figure 1: The strip of sky 9 degrees long from Fornax centre at the top left (NGC1399) out to NGC1291 at 
the bottom right. 
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niques. After the stacking and trimming, the size of 
each field was decreased to «1 deg. It is very im- 
portant for the automated detection of LSB and faint 
galaxies that the sky background is as flat as possible 
across each frame. Local background variations could 
lead to the detection of spurious objects in those areas 
where the sky appears brighter than its global value. 
In the opposite sense, real LSB objects would not 
be detected in regions where the local sky is fainter 



than the overall value ( Davies et al. 1994 ). To over- 
come this problem we have modelled the sky back- 
ground on each frame and then subtracted this from 
the data. The technique consists simply of creating a 
map of the background sky in which each pixel takes 
the median value of the pixels in the original frame in 
an nxn box surrounding the pixel in question. This 
map is then subtracted from the original data. The 
choice of box size is critical, as it will constrain the 
maximum size of objects that can be detected. Any 
object with a size approximately half or greater than 
the box size will be lost or at least severely degraded. 
The large scale size of the background variations over 
the frames allowed for a large smoothing box size. A 
340x340 pixel box was used. Considering the 2.3 
arc sec pixels of the CCD, this translates to objects 
with angular sizes of about 220 arc sec, very much 
larger than the objects we might expect to detect and 
our size selection criteria (section 2.2, below). Pre- 
vious work ( pavies fc Phillipps 1989| ) indicates that 
Fornax cluster LSB galaxies should have maximum 
diameters of about 50 arc sec. We will show below 
that all objects detected are much smaller than this 
upper size limit. This technique allowed us to flatten 
the sky to sa 0.5% which equates to a la fluctuation 
of «26 Rfi. 

Calibration was carried out using Landoldt stan- 
dard stars ( |Landolt 1992| ) and the standard Harris R 
filter used for the observations. The observing con- 
ditions were not ideal (some high cloud) and we es- 
timate a magnitude error of 0.2, which is sufficiently 
small for the comparison we intend to make. The 
sky brightness measured from the frames was ~20.5 
R/z, a value typical for the location, season and lu- 
nar age. Finally, astrometry was carried out using 
the ESA Guide Star Catalogue and is accurate to 
about one pixel, 2.3 arc sec. We used this astrometric 
calibration to measure distances from NGC1399 (the 
central cluster galaxy) and to obtain surface densities 



of galaxies in annuli around NGC1399. 



2.2 Image detection and object selection 

The image detection was carried out in two steps. 
First the frames were scanned for objects using 



the SExtractor package (Source Extractor, Bertin fc 
Arnouts 1996). All objects with 5 connected pixels 



brighter than 26 Rfi (the brightest 1 sigma isophote 
from all of the frames) were included in the catalogue. 
SExtractor also includes a neural network trained 
program that classifies each object with a "stellar- 
ity index" between (galaxy) and 1 (star). A total 
of 322,505 objects were detected in this way, giving 
an average of 23,200 objects per sq deg. The objects 
were detected by their size with the faintest having 
an isophotal magnitude of 22.8 Rfi. The SExtractor 
catalogue was then used as the input into our se- 
lection program this eliminated any flagged objects, 
such as ones close to the edges of the frame, near 
bright stars and composite objects. Also, all objects 
indexed 1 (definite star) were eliminated at this stage. 
The catalogue was thus restricted to non-flagged ob- 
jects with stellarity index up to 0.99 that met the 
selection criteria. With such large pixels it is diffi- 
cult to define a star, but the second stage selection 
described below shows how we finally selected only 
very extended objects. 

The final selection criteria are crucial because we 
wish to define a sample of galaxies that are predomi- 
nately nearby cluster members and not distant back- 
ground galaxies. Our first criteria was signal-to-noise 
ratio. To be sure we had real detections we set a min- 
imum signal-to-noise ratio of 10. This equates to a 
minimum area of ~ 28 pixels (148 sq arc sec). Using 
numerical modelling MDS showed that by selecting 
objects of LSB (in this case fj,g > 22.5) and with rela- 
tively large exponential scale sizes (a > 3" ) a cluster 
sample was predominately selected. This conclusion 
was also reached independently by [Schwartzenberg 
1996 and Phillipps et al. 1998 (we will return to the 



question of background contamination in section 3). 
We also set an upper central surface brightness limit 
of 23 Rfi so that we have a disjoint surface brightness 
distribution from the LSB sample (see section 3). So 
our final selection criteria for the VLSB sample were: 

1. Ho > 23 Rfi 
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2. a > 3.0 arc sec 



< \i >= fi - 2.5 log / + 5 log x - 0.75 



3. Isophotal area > 148 sq arc sec 

Plotting the isophotal magnitude of each object 
against the logarithm of its isophotal area we can 
draw lines of constant surface brightness and expo- 
nential scale size and select only those objects satis- 
fying the above selection criteria (fig. |2|). 

Using these criteria, 3462 objects were selected, an 
average of 251 objects per sq deg, less than 1% of the 
original number of detections. 

2.3 Photometry 

We have used the image parameters produced by 
SExtractor to determine the photometric properties 
of each object in the final catalogue. We will assume 
that the objects detected have exponential light pro- 
files ( pavies, 1989|) . 

The mean intensity of a detected object is 



Lis 
A is 



where Li so and Ai 
area respectively, 



< / >= 

, are the isophotal luminosity and 



L iso = 2irl a 2 f 



and 



Thus from the isophotal magnitude and area cal- 
culated by SExtractor we can calculate the central 
surface brightness /i Q (from above), the total magni- 
tude using 



m 



rrn 



+ 2.5 log/ 



and the scalelength using 



a 



10" 



It is worth commenting about the influence the 
large pixels may have had on our determination of 
the photometric parameters. We do not fit surface 
brightness profiles. Instead we derive the exponential 
parameters that correspond to the isophotal luminos- 
ity and area measured by SExtractor. Errors would 
arise if background pixels were included or excluded 
in the measured area because of the intrinsic fluctu- 
ations in the background. The background fluctua- 
tions (<t) are at about the same level as the limiting 
isophote so, the probability of including or excluding 
an extra pixel is at most 0.2 for Gaussian fluctuations 
(0.2™ for n pixels). Including or excluding one extra 
pixel to the area, and ignoring the small contribu- 
tion to the luminosity leads to an error in the surface 
brightness of only about 0.05 magnitudes even for the 
smallest galaxies detected. This will have little effect 
on our results. 



Ai 



irr. 



where the fraction of the total light contained within 
the limiting isophote is 



and 



/ = 1 



x = ^ 



Mo 



1.086 



fiL is the limiting surface brightness in magnitudes 
per sq arc sec and rz, is the limiting size at 
This leads to the following relationship between 
the measured mean surface brightness and the 
exponential central surface brightness. 



3 Results 

3.1 The three samples defined by their 
surface brightness and magnitude 

In what follows we will compare the clustering scale 
and the contribution to the total cluster light of three 
samples selected by their surface brightness. These 
are three progressively fainter central surface bright- 



ness samples. A bright galaxy sample from Jones & 



Jones 1980| , a LSB sample from MDS and our VLSB 
sample describe in this paper. We will also convert 
the VLSB sample surface brightnesses and magni- 
tudes into the B band using (B-R)=1.5, a value typi- 



cal of Fornax cluster dE galaxies ( Evans et al., 1990 ), 
so that we can compare this sample with the bright 
and LSB samples. 
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Figure 2: The isophotal magnitude - area plane (log) from one frame. Left: all detected objects. Right: 
Only selected objects. The two straight lines running diagonally mark the limits in surface brughtness. The 
lower is our central surface brightness limit of 23 R/x, and the upper is our threshold of 26 R/i. The curve 
identifies the 3" scalength limit and the vertical line the minimum area criteria. 



Jones & Jones list 64 Fornax cluster galaxies with 
m_e < 16. We have previously measured the sur- 
face brightness distribution of these galaxies ( pisney| 
et al., 1990| ) and we reproduce this in fig. 3. This 



sample illustrates the familiar Freeman result (Free- 
man, ] |970| ). The distribution of surface brightness is 
sharply peaked at a value of ~21.5 B/U with an rms 
width of w0.4 B/i. These galaxies will be used to de- 
fine the luminosity distribution as delineated by the 
bright galaxies. 

The Fornax cluster LSB sample, from the cata- 
logue of MDS, is a small part of a survey covering 
some 2400 sq deg of sky. The survey was carried 
out using UK Schmidt telescope photographic plates 
scanned by the Cambridge (U.K.) Automated Plate 
Measuring machine. All objectes were objectively se- 
lected to satisfy the selection criteria: 

1. /j >22.5 Bfi 

2. a>3.0 arc sec 

This sample is extensively discussed in MDS. The 
surface brightness distribution is disjoint from that of 



the bright galaxy sample with a peak at about 23.5 
(fig. 3). 

The VLSB sample has been discussed in section 2, 
and its surface brightness distribution is again shown 
in fig. ||. Note that fig. || should not be seen as 
a true distribution of surface brightness for galaxies 
over some 7 magnitudes. Each sample was selected 
in a quite different way, at different isophotal lev- 
els against different sky backgrounds. Each will be 
incomplete for galaxies with central surface bright- 
nesses close to the isophotal limit and the higher 
surface brightness galaxies in each sample will be se- 
lected against because at a given magnitude they will 
be smaller and so they will not satisfy the size crite- 
ria. We have made no attempt to correct the sam- 
ple for surface brightness incompletness because the 
true distribution of surface brightness is not known. 
Comprehensive discussions of how such corrections 



may be made are given in Disney & Phillipps 1987 



Davies 199C| and pVlcGaugh et al. 1995| . We will com- 
ment further on incompleteness in section 3.3. 

The three samples also show a clear trend in ap- 
parent magnitude with the lowest surface brightness 
galaxies having the faintest magnitudes. The mean 
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Figure 3: The B band surface brightness distributions of three different samples of galaxies in the Fornax 
cluster. The numbers in the VLSB sample have been scaled by a factor of 20 in order to fit on the same 
plot with the other two samples. 
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apparent magnitudes of the three samples are approx- 
imately 14.2, 18.1 and 19.8 respectively (M B =-17.1, 
-13.2 and -11.5 for a distance modulus of 31.3, Fer- 
guson 1990). 

3.2 The radial distribution of the three 
samples 

In this and the following section we will assume that 
all the galaxies in each sample are cluster members. 
We will justify this assumption at the end of this sec- 
tion. Using the three data sets described above we 
can measure the decrease in surface number density 
from the cluster centre, for galaxies of different sur- 
face brightnesses (fig. The number density data 
for the bright and LSB samples are taken from MDS. 
The data for the VLSB are binned over annulii of 0.25 
deg, with all radial distances measured from the cen- 
tral Fornax galaxy NGC1399. Fitting an exponential 
to the data we find scale sizes of 0.5 ± 0.1, 1.3 =t 0.4 
and 1.9 ±0.3 deg for the bright, LSB and VLSB sam- 
ples respectively (The central surface density of the 
VLSB sample is flat within the inner 1 deg and so 
we do not fit over this region. Including these points 
would give an even longer scale length). These corre- 
spond to approximately 0.16, 0.42 and 0.64 Mpc for 
a Fornax distance of 18.2 Mpc (Ferguson 1990). The 
LSB and VLSB samples are almost 3 and 4 times 
respectively more spatially extended than the bright 
galaxy sample. The result for the LSB sample (clus- 
tering scales some three times larger than the bright 
galaxies) has also been shown to be true using a much 
larger sample of galaxies, over a larger area of sky 
that includes other nearby galaxy groups (see MDS). 
It is also consistent with previous work on the spa- 
tial extent of more luminous LSB galaxies - they are 
less strongly clustered than the brighter galaxies (|Mo| 
st al. 1994). The exponential fit central galaxy num- 
ber densities are approximately 40 (HSB), 10 (LSB) 
and 800 (VLSB) per sq deg. This implies a sharp rise 
in the luminosity function at the faint end (slope of 
a — 2, see below). 






3.3 The contribution to the total luminos- 
ity of the cluster 

Having derived the parameters of the exponential fit 
to the number density distributions we can integrate 



Figure 4: The normalised surface density distribu- 
tions for the three samples. The top three are the 
distributions of the three different samples with pois- 
sonian errors. The bottom plot is the three samples 
overlayed, without errorbars for clarity. 
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to determine the total number of cluster galaxies ex- 
pected in each sample. This corresponds approx- 
imately to the total numbers we actually have in 
the bright and LSB samples. For the VLSB sam- 
ple the data only extend to large enough distances 
on one side of the cluster and so the integral takes 
account of the area not included. We can then use 
the mean luminosity of each sample to make a crude 
estimate of the ratio of total luminosities produced 
at each surface brightness. These are 1.0:0.1:0.5 for 
the bright:LSB:VLSB samples respectively. Adjust- 
ing the integral to account for the flattening of the 
VLSB number density within the central 1 deg makes 
little difference to the total luminosity (0.5 of the to- 
tal luminosity becomes 0.45). Integrating to 4° rather 
than infinity reduces the value from 0.5 to 0.3 of the 
total. Subtracting a background of 100 galaxies per 
sq deg (see below) alters the 0.5 to 0.3. If the clus- 
ter population ends at a radial distance of pa 4° and 
there is a background contamination of about 100 
galaxies per sq deg then the VLSB sample contributes 
about 0.1 to the total cluster luminosity. Over the in- 
ner 2° radius of the cluster the total luminosity ratio 
is 1.0:0.05:0.2 ignoring any background contribution. 
We will argue below that effectively all of the VLSB 
galaxies belong to the cluster and so we conclude that 
the VLSB galaxy population contributes a significant 
fraction of the total cluster light. 

It was our initial intention to use the range of lu- 
minosities in the VLSB sample to define the faint end 
of the luminosity function, but it is not possible to do 
this because of the way the sample has been selected, 
which itself was a consequence of trying to keep the 
background contamination to a minimum. The com- 
bination of the surface brightness and scale size limit 
forces a surface brightness size relation upon the data 
(see fig. H). At a given size we preferentially select 
faint galaxies because these are the ones that have 
low enough surface brightness. This would have had 
the effect of making the luminosity function appear 
steeper than it really is (this will affect other samples 
selected in the same way, |Schwartzenberg 1996 ). We 
cannot correct for this because we do not know the 
number of cluster galaxies that satisfy the individ- 
ual surface brightness limits, but fail to be selected 
because of their size. For example prinkwater et ak 



ness as the galaxies in the bright sample, but with 
much smaller sizes and fainter magnitudes. These 
were not included in the Jones fc Jones 1980| magni- 
tude limited sample. 

What we have done is to estimate a lower limit to 
the contribution made by the LSB and VLSB sam- 
ples to the luminosity function. We have done this 
by simply using the mean magnitude (standard de- 
viation is about 0.8 mag) and total numbers (inte- 
grated to infinity) of each sample. We have then nor- 
malised this to the total numbers of bright galaxies. 
In fig. |H| we have scaled our numbers so that the 
bright galaxy numbers correspond to the normalised 
luminosity functions given in [Smith et al. 1996| ( con- 
verted to the B band), pmith et al. 1996 show that 
the luminosity function of three more distant clus- 
ters steepens beyond Mb ~ —19.5 with a faint end 
slope of ~ —1.8. Our data (lower limits) is consistent 
with this steep slope and extends this result by some 
3.5 magnitudes to Mb ~ —12. It is consistent with 



1999(j]ave identified (using redshifts) Fornax cluster 
galaxies with approximately the same surface bright- 



a similar result obtained by [Phillipps et al. 1998| for 
the Virgo cluster. 

3.4 Contamination by non-cluster galax- 
ies 

The above result depends on there being little con- 
tamination of our sample by non-cluster galaxies. 
The radial surface density plots of course indicate 
a substantial number of cluster galaxies, but there is 
still the possibility of contamination of the sample by 
non-cluster galaxies. 

The bright sample is one where cluster membership 
has been assigned because of a measured redshift and 
so there is no ambiguity. Possible background con- 
tamination of the LSB sample has been extensively 
discussed in MDS. The conclusion was that at most 
20% of the galaxies could be in the background. 

To try and assess the background contamination of 
the VLSB sample we have analysed the surface den- 
sity of galaxies over the full 9° extent of the data. The 
data terminate with one field centred on the peculiar 
lenticular galaxy NGC1291 (fig.|l]). Our original mo- 
tivation for choosing this strip was that MDS had 
identified a number of relatively large LSB galaxies 
in the vicinity of NGC1291. A visual inspection of 
the field around NGC1291 on the CCD also clearly 
indicates a large number of LSB galaxies. 
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Figure 5: A comparison of our data with the luminosity functions derived for other clusters by Smith et al. 
1996. The filled circle indicates the bright galaxy sample. The LSB (open circle) and VLSB (star) samples 
are indicated as lower limits. Two Schecter LF are plotted (the fit to the Smith et al. data), with a = — 1.0 
(dotted) and with a = — 1.8 (solid). 
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In fig. |6] we show the normalised surface density 
of galaxies function of distance from the centre 
of Fornax extending to ~ 9 deg. The most strik- 
ing feature is the initial decline and then large in- 
crease in galaxy numbers as NGC1291 is approached. 
NGC1291 appears to be at the centre of a cluster of 
LSB low luminosity galaxies ! One can interpret fig. 
^ in two ways. Either the Fornax cluster population 
comes to an end at a radius of about 4 deg and there 
is a background surface density of order 100 galaxies 
per sq deg or the two extended LSB populations of 
Fornax and NGC1291 overlap at about 4 deg. 

Numerical simulations of the background contami- 
nation of the LSB sample by MDS predict about one 
background galaxy per sq deg with scale sizes greater 
than 3 arc sec. The same simulation but with the 
lower surface brightness cut-off of the VLSB sam- 
ple predicts zero background galaxies per sq deg (the 
simulation assumes a flat faint-end to the luminos- 
ity function ( Marke et al. 1994| ) and a Gaussian sur- 
face brightness distribution with a peak at Freeman's 
value ( Freeman 1970|) ). It is easy to show why this 
is so. A typical L* spiral galaxy would have to be 
at a redshift of z ~ 0.3 to have sufficient cosmolog- 
ical dimming to enter the VLSB sample. At z=0.3 
its scale size would be less than an arc second. To 
have an observed scale length of 3 arc sec it would 
need to be at z ~ 0.07 at which point it's surface 
brightness is cosmologically dimmed by only about 
0.5 magnitudes. The other option is that the field 
galaxy luminosity function is similar to that of the 
cluster ie. that beyond about 4 deg the galaxies in 
the VLSB sample are no longer associated with the 
cluster. We believe that there is observational evi- 
dence against this idea. 

Observations of the general field have been made 



by [Davies et al. 1994j and |Dalcanton et al. 1997| . Both 
of these surveys were optimised to detect LSB galax- 
ies in the field with similar surface brightnesses and 
scale sizes to the galaxies in the VLSB sample. The 



result from Davies et al. 1994 was an upper limit of 



10 per sq deg and from Dalcanton et al. 1997| 7 per 
sq deg. This is far fewer than the 100 per sq deg 
suggested by our data. So, our conclusion is that 
these are two overlapping extended haloes of dwarf 
LSB galaxies about Fornax and NGC1291. We be- 
lieve that the VLSB sample is essentially a Fornax 
cluster sample and that our three samples are not 



adversely affected by background contamination. 



3.5 NGC1291 

NGC1291 is not classified as a Fornax member. It 
is a peculiar lenticular barred galaxy in the field be- 
tween us and Fornax at a velocity of 839km s _1 ( |D"e| 
Vaucouleurs, 197lf ) . The mean velocity of the Fornax 



cluster is 1366km s 1 ( Ferguson, T990| ). In the same 
way as we fitted an exponential to the surface density 
of galaxies around Fornax we have derived the ex- 
ponential scale length of the surface density around 
NGC1291. The value is w 1.4 +/- 0.1 deg which 
equates to 0.22 Mpc at the distance of NGC1291 
(The Fornax scale length is 0.64 Mpc). This im- 
plies an extended, but not ridiculously large size for 
a halo of dwarf companions about an isolated galaxy 
( Zaritsky et al. 1993| ). We estimate the total lumi- 
nosity of the VLSB population to be about 0.7 that 
of NGC1291 (apparent B magnitide of 9.4). This is 
a marginally larger number than Fornax, but the ob- 
servations extend some 1.5 magnitudes further down 
the luminosity function and so we might expect about 
4 times more galaxies if the luminosity function car- 
ried on with a slope of -2.0, each with about 4 times 
less luminosity. This suggests that if you could iden- 
tify an ultra-LSB population you would conclude that 
there is more light in the extreme LSB constituents 
of the cluster than in the bright galaxies. Of course, 
at some point the luminosity function must turn over 
or terminate. 

It is interesting to calculate the approximate mass 
contained within the VLSB galaxy population sur- 
rounding NGC1291. pxwin &: Hatzidimitriou 1995 



find that typical (M/L) ratios for dwarf galaxies lie 
within the range 10-100, whilst that for the central, 
optical, part of a 'normal' galaxy is in the range 1-10 



(e.g. Kent 1987). Given our above estimate that the 
total amount of light emitted from the VLSB galax- 
ies is comparable to that from the central galaxy, 
over ten times as much mass may reside in the VLSB 
companions. This is comparable to the halo mass of 



2 x 10 12 M determined by |Zaritsky et al. 1993| for 



isolated spiral galaxies. Thus the agreement between 
the scale-length of the NGC1291 dwarf population 



and the halo size estimated by Zaritsky et al. 1993 



together with that between the halo mass and total 
mass of dwarfs, indicates that the 'dark halo' of some 
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Figure 6: The surface density distribution extending to 9° from the centre of the cluster. 
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spiral galaxies could be attributed to previously un- 
detected VLSB galaxies. 

The VLSB galaxies around NGC1291 may have 
a primordial origin, but another explanation may 
be found in the morphology of the bright galaxy; 
it is a peculiar galaxy. Its morphological classifi- 
cation ranges from "the remarkable barred lenticu- 
lar galaxy" ( Dc Vaucouleurs, 1975|) to a ringed SBa 
(NED) it has a dual personality because it also ap- 



pears to be known as NGCi^oy. It is quite possible 
that NGC1291 has relatively recently undergone a 
merger event and that some or all the VLSB galaxies 
are the debris from this collision. 



4 Discussion 

The existence of such large numbers of previously 
undetected galaxies is an important result and ones 
first thoughts are to other observations that might 
contradict it. One concern is that we do not seem 
to find very large numbers of low luminosity VLSB 
galaxies in the Local Group. The observed luminosity 
functi on of the Local Group is relatively flat over the 



range bf magnitudes sampled by our data (a ~ —1. 
van der Bergh 1992] ). Recently Blitz et al. 1999 have 
suggested that the Local Group mass function is re- 
ally quite steep. They have proposed that the high 
velocity clouds detected at 21cm are in fact the re- 
mains of the smaller dark haloes predicted by the nu- 
merical simulations. This remains to be confirmed. It 
may be that all groups have flat luminosity functions. 



Muriel et al. 1998 have compared the luminosity func- 
tion of galaxies in groups with those in clusters. They 
found that the luminosity function faint-end slope of 
the groups was much flatter (a ~ —1.0) than that of 
the clusters (a ~ —1.4). So the apparently flat lu- 
minosity function of the Local Group compared to a 
cluste r like Fornax may not be unusual or a concern. 



ated with the cluster low luminosity VLSB galaxies. 
Considering the selection effects and an assumed stel- 
lar population Theuns fc Warren 1997| estimate that 
these stars could account for 40% of the total clus- 
ter light a similar value to that derived by us for the 
VLSB galaxy population. 

There has been an on-going debate about the na- 
ture of the objects that give rise to the absorption fea- 



tures seen in the spectra of quasars (Lanzetta et al. 



1999). The issue is whether the Lyman limit and 



the damped Lya absorption features arise in nor- 
mal galaxies (in which case they have very large 
haloes) or in a previously unseen population of VLSB 
galaxies that are companions to the brighter galax- 
ies dPhillipps et al. 1993| , |Linder 1998]) . The Fornax 
VLSB population has a total cross-sectional area on 
the sky about three times greater than that of the 
brighter galaxies (equating isophotal sizes). If the 
VLSB galaxies have sufficiently high gas masses then 
one would expect the majority of quasar absorption 
features, seen in sight lines through the cluster, to 
occur in these small galaxies and not in extended 



haloes of the brighter galaxies. Recently Impey et al. 



1999 have compared absorption line features seen in 



With an observed central galaxy number density of 
about 450 per sq deg we predict a diffuse intra-cluster 
light of about 31 Bfi at the centre of the cluster. Such 
low surface brightness levels would be very difficult 
to detect against the relatively bright foreground sky 
(at best 23 Bfi). It is unlikely that a large popula- 
tion of low luminosity VLSB galaxies like this could 
be detected by their integrated light. Inter-cluster 
stars have previously been identified in Fornax by 



the spectra of quasars behind the Virgo cluster with 
known cluster galaxies. They point out the ambi- 
guity in trying to assign a given absorption feature 
with a particular galaxy. They also show that as you 
select Virgo galaxies of fainter magnitudes you are 
much more likely to conclude that a low luminosity 
galaxy is the absorber rather than a bright galaxy 
with a very large halo. [Impey et al. 199S also find 
that the correlation amplitude of the Ly Q absorbers 
is 4-5 times smaller than that of the brighter galax- 
ies, consistent with the absorption lines occuring in 
an extended VLSB galaxy population similar to that 
described here. Observations of Local Group dSph 
galaxies indicate low gas masses, « 10 7 M Q ( [Youngj 
fc Lo 1997 ), but because of the small size of the 



systems the column densities are not extremely low 



10 19 - 10 20 atoms 



cm 



This is typical of what one 



Theuns Sz Warren 1997j . These stars may be associ- 



might expect through the outer regions of a galac- 
tic disc. We conclude Ly a absorption lines could 
thus arise from an extended distribution of VLSB 
low luminosity galaxies in clusters or around isolated 
bright galaxies like NGC1291 rather than in extended 
gaseous halos. 

Others have measured steep cluster luminosity 
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functions before (e.g. Phillipps et al. 1998 , Tren- 
tham f997|, |Wilson et al. 1997| ). O ne concern about 
these determinations has been the possible contami- 
nation by background galaxies. The measured back- 
ground contamination was statistically removed from 
the data by assessing the numbers of background 
galaxies at each magnitude determined from fields 
away from the cluster. If this had been done in- 
correctly then a fit to the faint background num- 
ber counts (slope « 0.6) would lead erroneously to 
a steep (a rj —2.5) faint-end slope to the luminosity 
function. Our approach is different. We have specif- 
ically selected large LSB galaxies because our mod- 
els (MDS) indicated very low background contami- 
nation if we selected in this way. In addition because 
of the large extent of our data we have been able 
to show a decrease in galaxy number density with 
cluster radius. This strongly supports our contention 
that we have predominately selected a cluster sample 
and confirms the previous observations and interpre- 



tations of phillipps et al. 1998; , |Trentham 1997| and 



Wilson et al. 1997 



The clustering scale of the VLSB sample compared 
to the bright sample is similar to predictions of N- 
body hierarchical simulations of structure formation 
with "bias" ([White fc Frenk 1991| , |Frenk et al. 1996| , 



Kauffman et al. 1997). For example many features 
of the numerical simulations of a galaxy cluster by 



Frenk et al. 1996 are seen in the observational data. 



The galaxies show mass segregation with the most 
massive galaxies being more centrally clustered in re- 
lation to the smaller galaxies and to the dark matter 
(by about a factor of 3). Also the low mass slope of 
the simulated mass function is ~ —1.8. To be the 
dominant mass component of the cluster the VLSB 
galaxies require large values of (M/L). Previous work 
dMcGaugh fe de Blok 19981 , parignan fe~Fr eemaru 
1988j]i as suggested that LSB and dwarf galaxies are 
the most dark matter dominated galaxies of all. For 
example the mass to light ratios of LSB dwarf galax- 
ies in the Local Group can be large ((M/L)q=10- 
200, [Irwin fc Hatzidimitriou 1995|) . So, the VLSB 



galaxy population could account for all of the dark 
matter in clusters if the galaxy mass to light ratio 
is similar to the cluster mass to light ratio (of order 



100-150, David et al. 1995 ) or there are significantly 
more galaxies below our present detection limits. The 
VLSB galaxies may well be the most reliable tracers 



of the dark matter content of clusters and so we might 
expect the dominant dark matter of a cluster to re- 
side in its difficult-to-detect very LSB, low luminosity 
galaxy population. 

5 Conclusions 

We have identified a faint VLSB galaxy population 
associated with the nearby Fornax cluster. The 
VLSB galaxies appear to be associated with the 
brighter cluster galaxies, but their distribution is 
more spatially extended. Integrating the total num- 
bers we estimate that they contribute almost as much 
to the total luminosity of the cluster as the brighter 
galaxies and they have a larger cross-sectional area 
on the sky. If they have mass-to-light ratios as high 
as some more nearby dwarf galaxies they will domi- 
nate the mass of the cluster. Faint galaxies like these 
have previously been predicted by extrapolations of 
recently determined cluster luminosity functions and 
by recent numerical models of cluster formation. In 
the future it will be important to investigate the na- 
ture of this VLSB population further and to try and 
identify similar galaxies in other environments. 
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